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Abstract 

The Drinfeld-Sokolov reduction method has been used to associate with gl^ extensions 
of the matrix r-KdV system. Reductions of these systems to the fixed point sets of involutive 
Poisson maps, implementing reduction of gin to classical Lie algebras of type S, C, D, are 
here presented. Modifications corresponding, in the first place to factorisation of the Lax 
operator, and then to Wakimoto realisations of the current algebra components of the 
factorisation, are also described. 



1. Introduction 

We consider examples of constrained KP hierarchies having a Lax operator of the form 

L = e + z+{d + w)-^z_ with e = A'^d'' + uid''-'^ + --- + Ur, (1) 

where A is a constant diagonal matrix with having distinct, non-zero entries, Uq, ■ ■ ■ ,Ur 
are in glp, z+ G mat(p xs), Z- G mat(s xp) and w is in gig. Here mat(m x n) denotes the set 
of m X n complex matrices and for any vector space V, V stands for C°°{S^ , V) . PDO{m x n) 
denotes the set of pseudodifferential operators with coefficients in mat(m x 77.). We shall 
make use of the standard splitting PDO = PDO++PDO- of the space of pseudodifferential 
operators as a vector space direct sum of differential operators and integration operators. 
We also use the standard trace-form "res" on PDO given by resJ^ajQ* = a_i. 

We call the systems associated with Lax operators of the form given in (1) systems 
of extended G elf and- Dickey type. These systems are defined for any integers r,p > 1 and 
s > 0. For simplicity of language, we shall formulate our statements having in mind the 
generic case for which r > 1 and s > 0. Note however that all statements are also valid 
in the special cases for which either r = 1 or s = 0, even though some of them become 
trivial. The special cases for which r > 1 and s = reproduce the standard p x p matrix 
Gelfand-Dickey systems. The cases with r = 1 correspond to generalised AKNS systems. 

There have been several papers over the last few years devoted to systems of the above 
type [1-8]. It was shown in [9] how hierarchies with a Lax operator of the form in (1) can be 
obtained by the Drinfeld-Sokolov (DS) reduction method. Specifically, with the partition 

p times s times 

n = pr + s = r + - ■ + r + 1 -\ hi (2) 

is associated [11] a graded Heisenberg subalgebra of the loop algebra gin ® C [A, A~^], and 
a generalised KdV hierarchy having the Lax operator in (1) results from application of the 
DS reduction procedure [10] (see also [12,13,14]) with respect to a grade-one element from 
this Heisenberg subalgebra, if r > 1. In the r = 1 special case the DS reduction becomes 
trivial, but interesting results remain valid. 

The Lax operator usually studied in the literature for s > is obtained from (1) by 
choosing p — 1 and w — 0. In fact setting w = is not advantageous since this Dirac 
reduction of the phase space leads to non-local Poisson brackets. 

In the present paper we investigate the discrete symmetries given by involutive Poisson 
maps on the phase space of an extended Gelfand-Dickey system. Reduction to the fixed 
point set of such a map yields systems which arise from using the classical Lie algebras B, 
C, D in the DS approach. 

The further purpose of the paper is to study modifications of the above systems. In 
principle, modification arises via two possible mechanisms. The first can be viewed as an 
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application of the well-known factorisation approach of Kupershmidt- Wilson [15] (see also 
[16]). The phase space of the resulting modified system is a direct product whose factors 
carry linear Poisson structures, typically given by current algebras. The second is a novel 
construction which involves the so-called Wakimoto realisations of the current algebras, as 
was described in [17]. We will show that the two mechanisms are in fact closely related. 

We shall use the abbreviation "PB" for Poisson bracket and shall refer to the first and 
second Adler-Gelfand-Dickey PBs on PDO{p x p) [18,19] as the "AGD PBs". 



2. Extended Gelfand-Dickey hierarchies 

In this section we list the main elements of the theory of systems of extended Gelfand- 
Dickey type. Many of the results are described fully in [9] , whilst at the same time much of 
this theory is standard and goes back to the work of Adler, Gelfand-Dickey and Drinfeld- 
Sokolov, see [18], [19] and [10]. 

Let A4t)S^ be the space of quadruples (£, z-,w) that appear in (1), i.e., as a space 



■Mbs — [q^pY ^ niat(p X s) X mat(s x p) x gl^ 



(3) 



The functions on A^ds of interest are local functionals that have the form H = 
Jgi h{ui, . . . , Ur, z+, w) with h a difi'erential polynomial in the entries of the matrices in 
its arguments. There are two compatible PBs on A^ds, given by the following formulae for 
the respective hamiltonian vector fields: 
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Here the gradients are defined by 
d 



dt 



H{i+tSi, z±+tSz±, w+tSw) = Tr i^Si] + / tr 
=0 \oi J Js^ 
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(6) 



^ In the Hamiltonian reduction approach [9] A^ds represented a so-called DS gauge: we have 
kept this nomenclature here. 
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where Tr stands for J trres and ^ = ^2^=1 ^ ^1^- -^^^ ^ ^ PDO, we have the 
decomposition A = ^+ + A_ defined by the standard sphtting of PDO. 

The map tt : A^ds — PDO{p x p), which assigns the pseudodifferential operator L in 
(1) to the point {i, z+, Z-,w) in A4ds, is a Poisson map with respect to the PBs defined 
by the formulae (4) and (5) on A^ds and the first and second AGD PBs on PDO{p x p), 
respectively. It follows that M = 7r{M.Ds) — the set of Lax operators of the form (1) — is 
a Poisson subspace of PDO{p x p) with respect to the first and second AGD PBs. 

To define the commuting flows on M.-ds we proceed as follows. We first diagonalise the 
Lax operator L = ^{i, z^, Z-,w). That is to say we write 

L = gLg-' (7) 

for g an element of PDO{p x p) of the form g = Ip + 9kd~^ and we require that 

L be diagonal and g be off-diagonal, which determines them uniquely. For Q a constant, 
diagonal p x p matrix, let the functions H'^ be defined by 

HS{l,z±,w)^TT[LQ{^d)-'''), Hf{i,z±,w)^jTr (^P/^Q^ for j = 1, 2, . . . . (8) 

The set of functions for j = 0, 1, 2, . . . and Q arbitrary yields commuting Hamiltonians 
on A^DS- The corresponding hamiltonian vector fields are conveniently expressed in the 
following form: 



(^gQg-'L^/^) ,L 



X|,q(^+) = ^]+r,Qiz+) = res [gQg-'V/^z+id + w)'') 

X2q(^_) =X]+,^q(z_) = -res((a + «;)-iz_^Q^-iL^>) 

^iQi^) = ^]+r,QM =0' Vj = 0, 1, . . . . 

These commuting vector fields generate the fiows of the extended Gelfand-Dickey hierarchy. 
If r = p = 1 and s = 0, then the flows are trivial, and we henceforth exclude this case. 



3. Modifications of extended Gelfand-Dickey hierarchies 

We next apply a two-step factorisation procedure to the Lax operator L that leads to 
modifications of the flows in (9). By modification, we mean that there is a non-invertible 
Poisson map given in terms of a differential polynomial formula, from the Poisson space 
of the new (modified) variables to Air>s- The Hamiltonians of the modified flows are the 
pull-backs of the functions in (8). The flrst step of the factorisation procedure is 
rather well-known [4,5,6,9]. The second step was mentioned in passing in [9] but details 
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were omitted. Here we also explain the relationship of this second step to the Wakimoto 
realisations of the current algebra based on the general linear Lie algebra. 

Let us introduce the space = {glp) ^ ^ x glp^g and endow it with the current algebra 
PB on each of the components. The points of this space are denoted as (^i, . . . , 9r-i, Or) G 
0. For local functionals F^H on Q we thus have 



{F,H}{ei,...,er-i,er) = J2 [ 

i=i 



5F 6H 




(10) 



There is a Poisson map /j, from to Mbs described in [9]. It is important to note that 
is Poisson with respect to the second Poisson structure on A^ds given by (5), and not with 
respect to the first Poisson structure given by (4) . We shall not specify /i here, but we give 
the form of the composition $ = 7ro//:0^Mc PDO{p x p). 

Let us write the matrix 6r G glp+s ^^e form 

where a G glp, b G mat(p x s), c G mat(s x p), d & gig. Fix an integer k between and 
r — 1. Then $ is given by 

L = A(5+^i)A(5+^2) • ■ ■ A{d+e^)A[d+a-b{d+d)-^c]A{d+d^+-,) ■ ■■A{d+dr-i). (12) 

For K = there are no factors of the form {d + Oi) on the left, while for k, = r — 1 there 
are none on the right. Of course different choices of n correspond to different definitions 
of fi but all of them are related by invertible transformations. Hence all of the apparently 
different modifications for the different choices of re are equivalent. 

As the composition of two Poisson maps, $ is guaranteed to be a Poisson map with 
respect to the second AGD PB on PDO{p x p). A direct proof of the Poisson property of 
$ can be obtained using the following results. 

Lemma 1: The multiplication map : PDO x PDO — > PDO defines a Poisson map with 
respect to the second AGD PB on PDO. 

Lemma 2: \ 6 G glp} C PDO{pxp) is a Poisson subspace with respect to the second 

AGD PB, which on this subspace coincides with the current algebra PB appearing in (10) 
for i ^ r. 

Lemma 3: The map rj : glp^g PDO{p x p) defined by 

ry^^ ^^=d + a-b{d + d)-^c (13) 

is a Poisson map with respect to the current algebra PB on glp^g that occurs in (10) for 
i — r and the second AGD PB on PDO{p x p). 
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The first two lemmas seem to be part of the general knowledge in the field of integrable 
hierarchies. Lemma 3 was proved in [9] . The modified flows are defined on the phase space 
by pulling back the Hamiltonians in (8) by means of the map ji. 

For reasons explained in [9] (see also [1,2]), if s the factor 

K — d + a-hid + dy^c (14) 

entering the factorisation of L in (12) is called the "AKNS factor". Then the flows on 
can themselves be modified by factorising K as follows 

K ^d + a-h{d + d)-^c^{d + z?o) (Ip - l{d + ^?l + h)'^ 0) (15) 

for 

(t?o, 7) e glp X gig X S where S := mat(s x p) x mat(p x s). (16) 

We let the space glp x gl^ x S be endowed with the natural direct sum Poisson structure. 
That is if F, H are two local functionals on this space, we have 

6F rSHV^, r f5F5H _5H5F\ (17) 

The factorisation specified in (15) can be lifted to a mapping u : gl^ x gl^ x S —>■ glp+g 
whose equation is 

a = 1^0-^(3, d^'di + (3-f, b = i?o7 - 7^?i - 7/^7 + c = /3 (18) 
and direct calculation proves the following 

Proposition 4: If glp x gl^ x S and glp^g are endowed with the PB in (17) and with the 
current algebra PB, respectively, then the map v determined by (18) is a Poisson map. 

Define the space 0' = {glpY ^ x glp x gl^x S and endow it with the product PB given 
by the current algebra PB on {glpY ^ together with the PB in (17). Then ly gives rise to a 
Poisson map z/' : 0' ^ 0, which acts as u on glp x gl^ x S and as the identity on the {glpY 
factor. This map provides us with a modification of the system on = (glp)' ^ x glp^g. 
The resulting modified system is the same as the one engendered by the composite Poisson 
map n o I/' : Q' A^ds- 

We now explain that the map p defined by (18) can be used to generate a huge family 
of "realisations" of the current algebra PB based on gljn for any m. For this we simply 
repeat the construction for an arbitrary partition of m of the form m = mi + m2. This 



{F,H}{^o,^i,P,l) = 



5 



amounts to writing 9 e gl^ as ^ = ^ ^ J with a e glm^ etc, and expressing a, 6, c, o? by 
formula (18) in which we then insert the variables 

{'do,'di,P,'y) e gl^^xgl^^xSmi,m2 with cS^,,^^ := mat(m2Xmi)xmat(miXm2). (19) 

The PBs of these new variables defined similarly to (17) then imply the current algebra 
PB for the variable 6, that is we have a Poisson map like in Proposition 4. Repeating the 
construction iteratively for the current algebra factors, we can associate a Poisson map 

with any partition m = mi + 1112 + ■ ■ ■ + m/. The procedure gives that as a vector space 

Smi,m2,...,mi = ^'^ x ^'^ for 2d = u? - ml - 171^ mf, (21) 

and it carries the corresponding canonical PB. The precise formula of the map in (20) 
depends on the route whereby the final partition of m is reached through the iterative 
procedure. However, it has been proved in [17] (in a more general context) that the various 
Poisson maps that follow are all related by invertible Poisson maps. The map in (20) is 
known as a generalised (classical) Wakimoto realisation of the current algebra PB based 
on gljn- The standard Wakimoto realisation belongs to the partition of m for which all 
mi = 1. An explicit formula for the Wakimoto realisations was derived in [17] by different 
methods. Further results and background on Wakimoto realisations can also be found in 
[17] and references therein. 

We can use any of the Wakimoto realisations in (20) to modify any of the current 
algebra factors that appear in the factorisation of L in (12). This yields a large family of 
modifications of the extended Gelfand-Dickey hierarchies. 



4. Discrete reductions 



We now search for discrete symmetries of the extended Gelfand-Dickey systems. The 
compatible PBs on the phase space A^ds = Z-,w)} are given by 

{f,h}: = x^h{f) = 

for arbitrary local functionals F, H on A^ds^ where X^^^ are defined by (4), (5) for i = 1,2. 
Specifically, we look for symmetries given by some involutive map 

C7 :A1ds ^ A^DS, (T^^id, (23) 
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which leaves the PBs invariant, 



{Foa,Hoa}; = {F,H};oa, i = 1,2. 



(24) 



We take the following ansatz for a. Let m G GLp and let q G GLg, i.e., m and q are constant, 
invertible, respectively p x p and s x s matrices. Define the map (Jm,q '■ -Mds -Mds by 



/ £ \ / m£^m~'^ \ 



'm,q ■ 



Z- 



-mz^_q~^ 
qz\m~^ 
\ —qw^q~^ / 



(25) 



where is given by the standard adjoint operation on PDO{p x p), 

r r 

£t = (-lyA'-d'' + J2i-'^y~'9'-'ul for e = /Vd'' + Uid' 



(26) 



It is not hard to verify that (Trn,q satisfies (24) whenever it maps the phase space A^ds to 
itself, which is ensured by the condition 



mA^'m-^ = (-l)W. 



(27a) 



The involutivity of (Tm,q leads to the further conditions 

= e^m, = ±1, = Cg^, = ±1, with e^eg = -1. (27b) 

Notice that if = — 1 then p must be even and when = — 1 then s must be even. For 
any natural numbers a and h define the ax a and 2h x 26 matrices r]a and f226 by 



26 



h 2b 

^ei^26+l-i— ^ <^i,2b+l-i, 
i=b+l 



(28a) 



where the Cj are elementary matrices of appropriate size having a single nonzero entry 1 
at the ij position. Let denote an arbitrary a x a diagonal, invertible matrix subject to 



(28b) 



which means that is anti-symmetric under transpose with respect to the anti-diagonal. 
Using this notation, we have the following types of solutions for am,q- (The notion of a 
representative example is justified later in this section.) 

Type CI: r = 2p even, Vp, m is diagonal and q is arbitrary with Sg — —1, s — 21 even. 
Representative example: 0"A,n2r 



7 



Type C2: r = (2p + 1) odd, p = 2k even, A is such that rjpArjp = —A, m = ^pClp and q is 
arbitrary with eq = —1, s = 21 even. Representative example: crAQ2k,^2r 

Type B: r = (2p + 1) odd, p = 2k even, A is such that rjpArjp = —A, m = ^pijp and g is 
arbitrary with = +1, s = 21 + 1 odd. Representative example: CFAri2k,mi+i- 

Type D: r = (2p + 1) odd, p = 2k even, A is such that ripArfp = —A, m = ^pTjp and q is 
arbitrary with = +1, s = 2/ even. Representative example: (7Ari2k,V2r 

Note that the condition rjpArjp = —A, which is present except for type CI, requires A 
to have the form A = diag (Ai, . . . , A^, —A/-, . . . , — Ai), where p = 2k and A[ ^ i'^j 7^ 
for i ^ j since A^ must have distinct, non-zero entries. The Lie algebraic meaning of the 
notation referring to the various types will be explained below. 

Given an involutive symmetry a = am,qi one finds that cr : L i— > mL^m~^ for the 
Lax operator L in (1). It is not hard to see that this implies that the set of commuting 
Hamiltonians defined by eq. (8) admits a basis consisting of functions which are invariant 
or anti- invariant (that change sign) with respect to the action of a. On account of (24), 
ii H o a = H then the Hamiltonian vector fields are tangent to the fixed point set 
A4g,g C A^DS of a. Hence the flows of a "discrete-reduced hierarchy" may be defined by 
restricting the flows generated on A^ds by the cr-invariant Hamiltonians in (8) to the flxed 
point set A^gg. These flows are bihamiltonian with respect to the restricted Hamiltonians 
and a naturally induced bihamiltonian structure on Af^s- '^^e induced PBs on Al^s 
defined by restricting the original PBs of functions of cr-invariant linear combinations of the 
components of £, Zj^, Z-,w — which may be regarded as coordinates on A^£,s — to A^^g. 
The Lax operator of the discrete-reduced system belongs to 



For fixed p, r, s and a given symmetry type CI, C2, B or D the various possible choices 
of m and q defining am,q are equivalent from the point of view of the discrete reduction. 
In fact, the fixed point sets corresponding to two different choices are always related by a 
Poisson map of AIds given by 



with some constant matrices m G GLp and q G GLg. It is in this sense that the examples 
we gave for the symmetries of various types are representative examples. 

We mentioned that the extended Gelfand-Dickey system follows from an application 
of the DS reduction procedure to the Lie algebra gl^- As explained in particular cases in 
[20], the above discrete reductions are then induced by the reductions of gin to a simple 
complex Lie algebra Q of B, C or D type. This means that the discrete- reduced systems 



= n{M^s) = {LeM\L = mL^m-^ }. 



(29) 



{£, z+, Z-,w) t-^ {mim ^,mz+q ^,qz-m ^,qwq ^) 



(30) 
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are associated with graded semisimple elements of minimal positive grade from certain 
graded Heisenberg subalgebras of ^ (8) ([ [A, A"-*^] by means of DS reduction (see also [12,13] 
and the review in [14]). Since the graded Heisenberg subalgebras of Q <^ (L [A, A""*^] are 
classified [11] by the conjugacy classes [21] in the Weyl group W(^) of Q, we can label these 
generalised KdV hierarchies by the respective conjugacy classes in W(^). The conjugacy 
classes that occur here can be parametrised (as in [21,20]) by certain "signed partitions". 
The extended Gelfand-Dickey system itself belongs to the conjugacy class (r, . . . , r, 1, . . . , 1) 
of W{gln) given by the partition in (2). Using this notation, we find that the above discrete 
symmetries operate on the generalised KdV systems according to the following reduction 
rules: 

p times 21 times p times I times 

o-A,n2i '■ (2p, .r.,2p,l7'?77l) C W(5rZ2(pp+o) =^ (pT"^, iT*^^) C W{Cpp+i) 

2fe times ^SUimes^ fe times ^ I times ^ 

O-Afiafc.nji : (rv"^, l,-.-,l) C W(^Z2(fcr+0) =^ (rv"?^, 1, . . . , 1) C W(CA;r+0 
2k times 2 ^;+l^ime^ s times ^ I times ^ 

(^Ar,2k,V2i+i ■ (K^"^^. 1, . . . , 1 ) C W{gl2(kr+l) + l) =^ (^T^"^, 1, . . . , 1) C W{Bkr+l) 

2k times times^ times ^ I times ^ 

^Ar?2fc,r?2i : (^V'^^>1> •■•>!) C W(t/Z2(fcr+Z)) =^ {r^^^, 1, . . . , 1) C W (Dkr+l) 

(31) 

where Z > is arbitrary and r = 2p + 1 is odd. With the aid of case by case inspection, this 
result has been established in [20] for r > 1 and I = 0. Since the remaining cases can be 
treated in a similar way, we omit the proof (which simply amounts to diagram chasing). 

One may try to lift the discrete symmetries given in the above to analogous symmetries 
of the modified systems described in section 3. Considering the modified systems that 
correspond to L in (12), one needs to find a lifted transformation rule a- : — > for which 
/loa = ao fi, where : — > A^ds is the generalised Miura map. It is clear that such a local 
map exists if and only if the modification is symmetric in the sense that the same number 
of {d + 6) factors appears to the left and to the right of the special factor K in (12). This 
modification is available in the cases C2, B and D, for which r = 2p + 1 and we have 

L = A{d + ei)--- A{d + Op) A [d + a-h{d + d)-^c] A{d + Op+i) ■■■A{d + 02p) (32) 

by choosing k — p in (12). The transformation rule 9i ^ Of (0,1,..., 2p) is then not 
difficult to determine using the requirement that it must imply L ^ L" = mL'^m~^ for L 
in (32). Of course o" is a Poisson map, and the corresponding discrete-reduced hierarchy on 
the fixed point set 0^^ C provides a modification of the hierarchy on A^^g C A^ds- We 
leave it to the reader as an exercise to write down the explicit formula of a. 

For the discrete symmetry of type CI with / > 0, a factorised Lax operator of the 
symmetric form is available only after performing the second factorisation of K according 
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to (15). In this case r = 2p, and by choosing re = p — 1 in (12) (and renaming the variables) 
we indeed obtain the symmetric factorisation 

L = A{d + ei)--- A{d + Op) [Ip - 7(5 + ^ + /37)" V] A(a + ^p+i) ■ ■ ■ A{d + 02^). (33) 

The modified variables 9i {i — 1, . . . , 2p) and t?, /3, 7 now belong to the respective factors of 
the space 

e' = (glp)^^ X gl^i X mat (2/ x p) x mat(p x 21). (34) 

The lifted action of the discrete symmetry on ©' is easy to determine explicitly using that 
for L in (33) it must imply L >—>■ with a — crA,n2i • 

5. Concluding remarks 

We saw in section 4 that many KdV type hierarchies that are associated with certain 
conjugacy classes in the Weyl group W(^) for Q a classical simple Lie algebra by gener- 
alised DS reduction [12,13,14] are also obtained as discrete-reductions of extended matrix 
Gelfand-Dickey hierarchies. Note however that not all KdV type hierarchies based on a 
classical Lie algebra are discrete-reductions of hierarchies associated with gl^. For example, 
a pseudodifferential operator model of the KdV system associated with the primitive regular 
conjugacy class {p,p) in W(Z)2p) by generalised DS reduction is not known [20]. 

In the DS approach modifications of KdV type systems usually correspond to gauge 
transformations from certain "diagonal type gauges" parametrised by the modified variables 
to a "DS gauge" parametrised by the KdV fields. The map /i : G ^ TWds was obtained 
in [9] by this method. The modification p' : Q' ^ Q mentioned after Proposition 4 also 
permits interpretation as a gauge transformation in the DS approach. Moreover, the spe- 
cific factorisations of L in (32) and (33) that admit a local lifting of the relevant discrete 
symmetry have a clear interpretation. Namely, these modifications correspond to gauge sec- 
tions that are mapped to themselves by the original discrete-symmetry transformation that 
operates on the first order matrix differential operator variable used in the DS approach. 
More details on the way discrete symmetries occur in the DS framework can be found in 
[20]. 

Acknowledgement. The work of L.F. was supported by the Alexander von Humboldt 
Foundation. 
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